Excitation spectrum of the homogeneous spin liquid 



00 
On 
On 



O 

O 

oo 



> 

O 

00 
On 



c 

o 
o 



13 



R. Eder 

Institut fur Theoretische Physik, Universitat Wiirzburg, Am Hubland, 97074 Wiirzburg, Germany 

(February 1, 2008) 

We discuss the excitation spectrum of a disordered, isotropic and translationally invariant spin 
state in the 2D Heisenberg antiferromagnet. The starting point is the nearest-neighbor RVB state 
which plays the role of the vacuum of the theory, in a similar sense as the Neel state is the vacuum 
for antiferromagnetic spin wave theory. We discuss the elementary excitations of this state and show 
that these are not Fermionic spin-1/2 'spinons' but spin-1 excited dimers which must be modeled by 
bond Bosons. We derive an effective Hamiltonian describing the excited dimers which is formally 
analogous to spin wave theory. Condensation of the bond-Bosons at zero temperature into the state 
with momentum (n, n) is shown to be equivalent to antiferromagnetic ordering. The latter is a key 
ingredient for a microscopic interpretation of Zhang's SO (5) theory of cuprate superconductivity. 
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I. INTRODUCTION 

High-temperature superconductivity occurs in a state 
which is is frequently referred to as an 'RVB spin liq- 
uid'. This state has no magnetic order, but strong short 
range antiferromagnetic correlations. Undoubtedly the 
strong repulsion between electrons, which turns the sys- 
tem into a charge-transfer insulator at half-filling persists 
in the doped case, so that a description in terms of free- 
electron like Slater-determinants is not really adequate. 
The problem then is, how to describe such a state theo- 
retically. Despite its frequently being referred to in the 
literature, the RVB spin liquid is a rather elusive con- 
cept. For example the precise nature of its ground state 
and low lying elementary excitations is not known to any 
degree of certainty. In the following we want to address 
this problem by studying a disordered state for the 2D 
Heisenberg antiferromagnet 



H 



on an 2D square lattice. Here Si denotes a spin 1/2 op- 
erator on site i. One might expect that this is a kind of 
stepping stone also for the doped case, in that the ele- 
mentary excitations of the undoped spin liquid persist to 
some degree also for finite doping. 

Perhaps the best-defined RVB spin liquid is the nearest 
neighbor RVB state jl]-fl| - at least this wave function can 
be written down in compact form. We define the singlet 
generation operator on the bond i,j 
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where ^ =0^0^ 



are the constrained Fermion op- 
erators, which do not allow the creation of a second elec- 
tron on an already singly occupied site. Introducing the 
operator 
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where e.g. i + x denotes the nearest neighbor of site i 
in ^-direction, the nearest neighbor RVB state on a 2D 
square lattice with 27V sites can be written as 
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where n is a normalization factor. It corresponds to a 
superposition of all states which can be obtained by cov- 
ering the plane compactly with nearest neighbor singlets, 
all with equal phase. Covering the plane with singlets 
is equivalent to covering it with dimers, a well-known 
problem from statistical mechanics || . We can therefore 
rewrite the state as 
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where a denotes a dimer covering of the lattice and \?p a } 
the state obtained by putting singlets onto the dimers of 
a. 

In the following, we want to examine the problem of the 
possible elementary excitations of such a singlet back- 
ground, and set up an effective Hamiltonian governing 
their dynamics. The idea of 'expanding' around a suit- 
ably chosen vacuum state is realized in simplest form in 
linear spin wave theory. The general line of thought here 
is quite analogous to linear spin wave theory, with the 
sole exception that the role of the vacuum (which deter- 
mines the symmetries of the ground state) is played by 
the 'singlet soup' (|J) instead of the Neel state. A similar 
approach has previously been applied to dimerized pla- 
nar Heisenberg-type models [||J^] to spin ladders |,|, to 
strongly coupled Heisenberg planes |l(]|ll| , and to Spin- 
Peierls-like spin chains jl2|,[l3]]. An example where the 
fluctuations are Fermionic rather than Bosonic in nature 
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is provided by the Kondo lattice JT40. The main differ- 
ence as compared to the present work is that in all of 
these works a rather unique and simple dimer covering 
of the system is given by the topology or the form of 
the Hamiltonian - the complications that arise from the 
use of a disordered 'singlet soup' such as (||) then can be 
avoided. 

While the technical complications arising from the use 
of a dimer basis are considerable, this approach also has 
some major advantages: in a site basis it is virtually 
impossible to even write down a disordered spin state, 
because one has to deal with the spin degeneracy on 
each single site. The calculation only becomes feasible 
if this site-degeneracy is lifted, for example by assuming 
strict Neel order as in spin wave theory. On the other 
hand, the degeneracy is automatically taken care of in 
the dimer basis, because two interacting sites do have a 
unique ground state. A further considerable advantage of 
the dimer basis is, that it is easily enlarged by hole pairs 
on nearest neighbors, so as to describe a superconducting 
state. Indeed, as will be seen below, the present descrip- 
tion of the antiferromagnetic phase most naturally can 
be generalized to comprise also a superconducting phase, 
thereby providing a very simple microscopical picture for 
the SO (5) rotations which smoothly connect antiferro- 
magnetic and superconducting phase in Zhang's theory 
|1| of cuprate superconductors. 



II. ELEMENTARY EXCITATIONS OF A 
'SINGLET SOUP' 

The nearest neighbor-RVB state ([l|) has the symme- 
try properties expected for a homogeneous spin liquid: 
it is isotropic, translationally invariant, is an exact spin 
singlet and has no magnetic order. On the other hand, 
just as the Neel state, it is not an eigenstate of H. If we 
take one singlet configuration \ip a ) and act with the ex- 
change term on a bond connecting two different singlets 
(see Figure [|a) we can create a state which no longer can 
be represented as a superposition of only nearest neigh- 
bor singlets. Such a state therefore represents a kind of 
fluctuation and as a first step we need to understand the 
character of these fluctuations. It might appear |l6| that 
the energetically most favorable fluctuation is the state 
shown in Figure [j]b: two nearest neighbor singlets are 
transformed into a configuration with only one nearest 
neighbor singlet and a second singlet connecting more 
distant sites. Nominally the energy increases by only 
3 J/4 in this transition, because the only change is the 
loss of one nearest neighbor singlet. Because singlet and 
triplet are degenerate for sites which are not connected 
by an exchange bond, we might as well consider the two 
spins connected by the 'long singlet' as being unpaired 
(see Figure |l]c) . The transition from Figure [l] a^b could 
thus be viewed as pair creation of two unpaired spins. 



Next, by acting with the exchange term on a bond which 
connects a dangling spin to another singlet (see Figure 
|l|c) , we can recouple the spins so as to form a new singlet 
and leave one of the formerly paired spins unpaired (see 
Figure [lid) . This process corresponds to a propagation 
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FIG. 1. An incorrect picture for fluctuations in a singlet 
background. 

of the dangling spin. We would thus arrive at the con- 
clusion that the fluctuations out of the nearest neighbor 
RVB state are unpaired spins, which carry a spin of 1/2 
and consequently must obey Fermi statistics. Clearly, 
these excitations should be identified with the ominous 
'spinons'. 

Further reasoning shows, however, that the line of 
thought leading to the introduction of the 'spinons' 
misses a small but crucial detail. The first reason is that 
the state in Figure is not orthogonal to the vacuum, 



Piefactor: -1/2 
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FIG. 2. Placing singlets in either of the two ways shown on 
the left produces the spin configuration on the right with the 
indicated prefactors. This gives a contribution of —1/4 to the 
overlap, an equal contribution comes from the spin reversed 
configuration. 

and thus cannot represent a true fluctuation. More pre- 
cisely, if we introduce (see Figures [l] and g)) 
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it is straightforward to see (Figure |^) that (a\b}=— 5, 
in other words: after the 'transition' Figure |l]a— >b we re- 
main in the original state, Figure [lja, with a probability of 
25%. The problem of non-orthogonality is not restricted 
to the first step in Figure [1} the states Fig. 0c and 0d 
have an overlap of 1/2, and this generalizes to any two 
states which differ by one hopping process of a 'spinon'. 
The non-orthogonality problem thus is omnipresent and 



severe. 



Let us therefore return to the first step, Figure |l] a^b, 
and consider how we can remedy the problem. The most 
natural way to proceed is to form the orthogonal com- 
plement 

\b') = \b) - (a\b) \a), 

so as to see 'what is really new' in the state \b). A 
straightforward computation shows that after normaliza- 
tion to unity the orthogonal complement is 
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Here we have introduced the operators |17|,| 
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which create the three components of the triplet on the 
bond (i, j). We arrive at the conclusion that the true fluc- 
tuation out of the nearest neighbor singlet background is 
not the creation of two Fermionic 'spinons', but rather 
the creation of two excited dimers, which carry a spin of 
1 and consequently should be modeled by bond-Bosons 
The further evolution of the created triplets then 
is quite obvious (see Figure ||) (but completely different 
from that of the 'spinons'): by exchange along bonds con- 
necting the triplets with neighboring singlets the triplets 
can de-excite while simultaneously the singlet turns into 
a triplet - this process, which is very much reminiscent of 
the propagation of a Frenkel-type exciton corresponds to 
the propagation of the excited dimer. Note that unlike 
the 'spinon' states in Figure [l|, all different states in Fig- 
ure [| are mutually rigorously orthogonal. As a matter of 
fact there are problems with non-orthogonalities also for 
the triplet states - these are 'inherited' from the original 
nearest neighbor RVB state. They will be discussed in de- 
tail below and be shown to be much less severe than those 
for the 'spinon' states. Their main effect is to replace the 
simple excited dimer by a more delocalized object, which 



resonates between different orientations within a limited 
spatial region. 
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FIG. 3. A more correct picture for fluctuations in a singlet 
background. 

To be more precise, we now discuss the action of the 
Heisenberg exchange on all possible configurations of 
nearest neighbor singlets and triplets. Consider two near- 
est neighbor bonds and and assume that 
they are connected by a single bond (i, i'). Denoting 
the Heisenberg exchange along the latter bond by hi.i> 
we have: 
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These equations show that if we start out from states con- 
taining nearest neighbor singlets or triplets on the l.h.s., 
the exchange term only produces states which again con- 
sist of nearest neighbor singlets or triplets on the r.h.s 
(had the two bonds been connected by exchange along 
another bond, {hi') or (j,j r ), we would have ob- 

tained the same equations with the sole difference that in 
some cases the prefactors change their sign) . This proves 
rigorously that acting with an arbitrarily high power of 
the Hamiltonian onto the nearest neighbor RVB state 
produces only states which can be built up from nearest 
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neighbor singlets or triplets. 

This 'theorem' in fact holds true in a more general sense: 
on a single dimer to, the 4 states and do form a 
complete basis [fjj. Thus, if we use a fixed dimer cov- 
ering a, the set of states obtained by placing singlets 
and triplets on the dimers of a form a complete basis 
of the Hilbert space. Adding up such states obtained 
from all possible dimer coverings then clearly produces 
a highly overcomplete basis of the Hilbert space, which 
therefore must automatically include states with singlets 
of arbitrary length. It follows that all states with longer- 
range singlets also can be represented as superpositions 
of states which are composed exclusively from nearest- 
neighbor singlets and triplets. These states are therefore 
redundant, and if we formulate a self-consistent theory in 
terms of nearest-neighbor singlets and triplets, we have 
automatically included these longer-ranged singlets. The 
fact that we are using a nearest-neighbor RVB state as 
the starting point for constructing our theory therefore 
means no loss of generality and in particular does by 
no means imply that we are considering only states with 
only very short-ranged antiferromagnetic correlations. In 
fact, it will be shown below that one can construct even 
states with infinite-range antiferromagnetic order by us- 
ing exclusively nearest neighbor singlets and triplets. 
The preceding considerations suggest that we should 
model the excitation spectrum of the nearest neighbor 
singlet vacuum by Bosonic excitations, which approxi- 
mately correspond to excited dimers. Assuming that the 
bonds in the system have been labeled in some way, we 
denote the triplet operator on bond i by tj. Then, we 
introduce the following basis states 
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where n(iia±, iicti-, ■ ■ ■ i m oim) is a normalization factor. 
They describe a certain number (to) of triplets which are 
'immersed into the singlet soup'. Thereby the singlets fill 
the space in between the triplets compactly in all possible 
ways. All states which can be generated by pair creation 
and propagation of triplet bonds (such as the ones shown 
in Figure ||) can be represented in this way. We next con- 
sider the triplets as Boson-like elementary excitations of 
the singlet-background, in precisely the same way as mis- 
aligned spins are considered as Bosonic excitations in a 
'Neel background' in antiferromagnetic spin wave theory. 
Re-interpreting the state 



H = jY j r\-T i +Y J ^ijr\-T] + H.c. 
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where we have grouped the three triplet components into 
a 3- vector r so as to stress manifest rotation invariance. 
The first term in (^) corresponds to the energy of forma- 
tion of the triplets, the second term describes pair cre- 
ation processes as in Figure ||a — > ||b, and the third term 
accounts for the propagation of the triplets, see Figure 
||c — > ||d. The matrix elements and Ay should be ob- 
tained by computing matrix elements of the Heisenberg 
Hamiltonian H between the corresponding states (|^) . Of 
course one thereby has to assume that for example the 
matrix element for a triplet jumping from bond to to 
bond n does not depend significantly on the positions of 
the other triplets - otherwise a description in terms of 
a single-particle like Hamiltonian would not be feasible. 
As is the case in spin wave theory, the r-Bosons have to 
obey a hard-core constraint, and in fact presence of one 
Boson prohibits the presence of another Boson not only 
on the same bond, but also on all bonds which share a 
site with the original one. 

In the following, we will first study the problem of a 
single excited dimer in the singlet background, in other 
words we want to compute the 'bare' Boson dispersion 
e(k) in (^|). As our basis states we consequently choose 
(suppressing the x y or z spin-index of the triplet): 
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where a—x, y denotes the direction of the bond in real 
space, and ri\ a normalization factor. In this state one 
triplet is put onto the bond (i, i+a) and the remainder of 
the lattice is covered compactly by singlets in all possible 
ways. Next, we introduce the Fourier transforms 



Ife, a) — 
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In the Hilbert space of bond-Bosons, this state would be 
denoted by a |0). The procedure to be followed then is 
like this: in a first step we compute the 2x2 overlap ma- 
trix N(k)—(T k a ffc a i) and diagonalize it. Denoting the 
resulting eigenvectors and eigenvalues by e„ and A„, the 
states 



where the r- represent Boson operators, we may ex- 
pect to describe the dynamics of these Bosons by a 
Hamiltonian of the form flial 
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form an orthonormal basis set and hence can serve as ef- 
fective single particle orbitals with momentum k. Since 
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the Boson operators which correspond to the original 
triplets obey Q [f^f^] = N a ,p, the operators Tk, v 
obey the canonical commutation relations for Boson op- 
erators: [i)e,i/)4i/'] = &v,v' ■ They describe a triplet- 
like excitation which oscillates between x and y-directed 
bonds within a certain spatial region whose extent is de- 
termined by the range of the real space overlap integrals 
(i,a\j,a'). 

Next, we set up the 2x2 Hamilton matrix 



H(k)=( 



which in turn requires knowledge 



of the real space matrix elements (i, a\H\j, a'}. Diago- 
nalizing H(k) — EoN(k), where Eq denotes the expecta- 
tion value of H in the 'background' state (|J), we obtain 
the desired dispersion relation of a single triplet Boson. 
The pair creation matrix element is obtained in an anal- 
ogous way 

This procedure in fact is neither new nor unconventional: 
a completely analogous construction is performed e.g. in 
the construction of the t — J model |^(| , which describes 
the dynamics of the (non-orthogonal) Zhang- Rice singlets 
on the different plaquettes of the CuC>2 plane. The only 
difference is that here we have two different nonorthog- 
onal objects (Bosons on bonds in x and y-direction) per 
unit cell, whereas it was only a single Zhang-Rice sin- 
glet/unit cell in the case of the CuC>2 plane. Apart from 
that the construction is precisely the same. 
In the next three sections we will calculate the disper- 
sion relation, the pair creation matrix element, and dis- 
cuss how these matrix elements depend on the density of 
triplets. Readers which are not interested in these more 
technical parts are advised to proceed to section VI. 



III. PROPAGATION OF A SINGLE TRIPLET 

To carry out our program we need to compute the real- 
space matrix elements (i, a\j, (3) and (i, a\H\ j, (3). In do- 
ing so the concept of a loop covering of the plane Q is of 
crucial importance. For two dimer coverings a and b the 
loop covering c = a + b is obtained by drawing a and b 'on 
top of each other' (see for example Figure 1 in Ref. 0). 
This produces a covering of the plane by closed loops 
u, each of even length 2L(u) (note that in the following 
we always measure the length of a loop L(u) 'in units of 
dimers'). Let us now consider each loop u as an isolated 
ID ring with 2L(u) sites. We assume that the sites along 
the ring are labeled such that the dimer covering a corre- 
sponds to the state \a) — Y\i s\ i+1 \0) - the dimer covering 

b then must correspond to \b) = ]J i s\ +1 i+2 \ty- Expand- 
ing the products we get 2 L ^ different spin states from 
each covering, and there are precisely 2 spin states which 
show up in both \a) and \b) namely the two possible Neel 
states. We thus have (a\b) = 2/2 L( - u \ The same holds 
true for any other loop, whence, using X)«ea+6 L(u) = N, 
we find the scalar product of the two singlet distributions 



? P(a+b)-Af 



where P(c) is the total number of loops in the loop cov- 
ering c. 

Let us now assume that the singlets on the bond m in 
\ip a ) and on the bond n in \ipb) have been replaced by a 
z-like triplet (due to the explicit rotational invariance of 
the 'singlet soup' the result for an a:-like or y-like Boson 
will be precisely the same - we are choosing the z-like 
component because the ambiguous states in this case are 
again the ones with Neel order along the loop). Then, a 
necessary condition for the scalar product to be different 
from zero is that there is a single loop uq in the resulting 
loop covering a + b which passes through both bonds m 
and n. The reason is that the time-reversal parity of the 
triplet is negative whereas that of the singlet is positive. 
A necessary condition for a loop u to give a nonvanish- 
ing overlap is that the total time-reversal parities 'along 
the loop' are equal for both states \tp a ) and \ipb)- This, 
however, is only possible if the triplets in \i/; a ) and \ipb) 
are within the same loop. Each loop in a + b therefore 
must contain either no triplet or both of them. 
We can now split up the entire overlap integral (ip a \ipb) 
into components which differ by the length and topol- 
ogy of the loop uq which passes through both triplets. 
The absolute numerical value of the overlap from this 
particular loop is identical to the case of pure singlet 
covering. The only change may be an extra minus sign, 
which originates because the singlets do have an orien- 
tation, whereas the triplets do not. We thus can rewrite 
the overlap as 

(Tp a \i> b )=J2 2 1 - L( " o) (-1)<**°) X M, 

2 -(N-L(u)) 



X(u) 



ni 



£ 2 P(a +b )-l Aa+M _ (13) 



Here A c u is 1 if the loop covering c contains the loop 
u and zero otherwise. We also note that %(1)=1, which 
fixes the normalization factor n\. With the exception of 
the x( u ) all parts in ( [l3|) can be computed analytically. 
x(u) may be viewed as the norm of a nearest- neighbor 
RVB state which covers only the exterior of the loop u, 
divided by the norm n\ of the state which covers the 
exterior of a single bond. If we assume that the norm 
increases exponentially with the number of sites in the 
system, n i=s e aN , with a > 0, one would estimate that 
x(u) ~ e~ Q( - L( -"' l ~ 1 ). This suggests that x( u ) is a quite 
rapidly decreasing function of L(u). In the present work 
numerical values for x( u ) were obtained by exact calcu- 
lations on finite clusters (see the Appendix). The x( u ) 
thereby indeed turned out to decay rapidly with L(u), so 
only contributions with L(u)<2 were kept in the present 
calculation. It should be noted that the computation 
of the x( u ) is n0 fundamental obstacle to the present 
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scheme: it may well be possible to obtain essentially 
exact values for these parameters by using Monte-Carlo 
techniques on large lattices. Figure then shows the 





(a) 



FIG. 4. Pairs of bonds which are connected by L=2 loops 
(a) and L=3 loops (b). Bond i is kept fixed (dashed line), 
bond j (full line) is labeled by (-1) ct(u) . 

pairs of bonds (i, a) and (j,f3) which can be connected 
by loops of length 2 and 3 as well as the corresponding 
signs (— l) CT ( Uo ). In this way we find the overlap matrix 



with 7i(fe) = l and 



72 (ft) = 







^ 4sin(^)sin(^ 



4sin(^)sin(- 







(14) 



We proceed to a calculation of the matrix elements of the 
Hamiltonian. We first recall that for the nearest neighbor 
RVB state the expectation value of H between two dimcr 
coverings \ip a ) and \ipb) can be decomposed into contri- 
butions from each individual loop in the loop covering 
a + b§: 

u£a-\-b 

E(u) = e s [(2 - <$l(„),i ) L(u) + n b (u) ], 

where 7i&(ti) is the number of nearest-neighbor bonds in 
u which bridge the loop (see Figure ||), and e s =— 3J/4 
the exchange energy/singlet. E(u)/e s is the number of 
nearest-neighbor bonds that can be formed from the sites 
covered by u 0] . This formula implies that there is no 



(h - - - -£) 



contribution from bonds connecting different loops. The 
reason is that the exchange along a bond connecting dif- 
ferent singlets can only lead to the pair creation of two 
triplets, see the first equation (0). In order to maintain 
time reversal symmetry along each loop, it is then neces- 
sary that both triplets belong to the same loop - which 
is not possible if the bond in question connects different 
loops. 

Let us now again assume that the bond to in \ip a ) and the 
bond n in \ipb) are occupied by a triplet, and consider the 
'connected matrix element' of H between the two result- 
ing states: (ipb\H\ip a ) — Eo{ipb\ip a ). First, let us assume 
that we act with the exchange along a bond connecting to 
and a neighboring bond to'; the triplet can either prop- 
agate from to to to', or decay into two triplets on both 
to and to' (see the second Equation @)). Neglecting the 
second possibility we obtain a nonvanishing contribution 
to the matrix element of H only if there is a single loop 
Mo € a + b which covers both, to' and n. Alternatively, 
if we act on a bond which connects n and a neighboring 
bond n' , the triplet jumps from n to n and we get a non- 
vanishing contribution only if one single loop uq G a + b 
covers n! and to. If, on the other hand, we act with the 
exchange along a bond which does not touch either of 
the triplet bonds to or n, both triplets will stay where 
they are and we get a nonvanishing matrix element only 
if both, to and n, are covered by a single loop uq a + b. 
The same holds true if we act with the exchange along 
the bonds to and n themselves. 

We first consider the case that to and n belong to two dif- 
ferent loops, Mo and u' . In the simplest case both 'loops' 
consist only of a single bond, i.e. Uq consists of the single 
bond to and u' Q only of n. Since the two triplets belong to 
different loops, the overlap (ipb\i/) a ) is zero. Moreover, the 
matrix elements of the exchange along any bond which 
does not connect to and n vanishes - the calculation thus 
becomes very easy. The matrix element for the triplet 
hopping from to to n is ±J/4, where the signs for dif- 
ferent relative positions of the two bonds are shown in 
Figure |(| To 'embed' the hopping process 
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FIG. 5. 'Bridging bonds' (dashed lines) in the L — 3 loop 
(a) and in an L = 4 loop (b). 



FIG. 6. The sign of the hopping integral from the dashed 
bond m (dashed) to the indicated bond n. 
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into the singlet background, we need to renormalize this 
matrix element by 

9-JV 

Vnm = V 2 p ( Q+b ) A 0+6>ro A a+6>n , (15) 

ni ^ 

which is again estimated from cluster calculations. The 
first contribution to the Hamilton matrix then is: 

e (D (fe ) = il h (k), (16) 

where the elements of the 2x2 matrix ti(k) are: 

ti,xx(k) = 4 cos(k y ) — 2 cos(2fc.,,) — 4 cos(fcj.) cos(fcj,), 
ti :V y(k) — 4 cos(k x ) — 2 cos(2fcj,) — 4 cos^) cos(k y ), 

^ y (fc) = 4sin(^)sin(^)+4sin(^)sin(^). (17) 

To keep things simple we have moreover replaced the dif- 
ferent r) nm by the average value r\. 
Next, we consider the case that one of the loops, e.g. u' , 
has a length > 2. In other words, we a considering a 
process like the one shown in Figure ^ the triplet jumps 
from bond m to bond m', and the triplet on m! 




(a) (b) (c) 

FIG. 7. By application of H a triplet can hop from bond 
ra — > m! (see (a)— ^(b)) and then overlap with the triplet in 
the final state on bond n along the indicated L = 2-loop (see 
(b)^(c)). This process gives a nonvanishing hopping matrix 
element from bond m — > n. 

overlaps with the triplet on bond n along a loop (in this 
case of length 2). There is also an analogous process, 
where the triplet jumps from n to m! and m! and m 
overlap by a loop. The respective matrix elements can 
be factorized into the matrix element for the hopping of 
the triplet times the overlap along the loop. This means 
that the matrix element can be written as ^^(to, n ) 
where 

t 2 (m,n) = ^[ti(m,TO')72(m',n) + ti (n, m')j 2 ('m' 1 m)]; 

ra' 

here t\(m, m') and 72 {m! , n) are the real-space versions 
of the matrices ( |l7| ) and (pd|). Fourier transformation 
gives 

c (2)( fe ) = ^ l( fc) 72( fc) +72 (fc)t l( fc)) (is) 



where it has to be kept in mind that ii(fc) and 72(&) are 
symmetric 2x2 matrices. The 'embedding factor' is 

o-iV 

m = V 2 p ^A a+b . m A a+b . n ,A a+b . n (19) 

where (m, n', n) are like in Figure |^. Actually there are 
two inequivalent relative orientations of a single bond m 
and an L = 2 loop - the two respective values of 773 do 
not differ strongly and for simplicity we use the average 
of the two values for the two configurations. Processes 
involving even longer loops could be treated in an anal- 
ogous way, but we neglect these. 

We proceed to the second case, i.e. we assume that m 
and n are covered by a single loop uq. As was the case 
for the overlap integrals, the matrix elements then can be 
split up into contributions differing by the loop Uq which 
covers both triplet bonds. Once this loop is fixed, we can 
divide all bonds in the plane into three distinct classes. 
First, bonds belonging to any loop other than uq are not 
affected at all by the presence of the triplets and give 
the same contribution as in the pure singlet state. This 
becomes 

2 -(A f -L(«o)) 
2l-L(u ) (_ -|\o"(«o) 

Til 

Yy E ^(")]2 p(a+b) - i A Q+b , U0 . (20) 

a,b u e a + b 

This is an energy of order of the system size N; in the end 
it must be canceled, up to terms of order TV , by a corre- 
sponding contribution in —Eo(tf> a \ipb). This cancellation 
is the analogue of the familiar linkcd-cluster theorem of 
many-body physics. 

Second, we consider those bonds in uq which are not 
bridging bonds. They will be covered by either a singlet 
or a triplet in either a or b, whence these bonds together 
give the contribution 

[E(uq) + (2 - h(u a ),i )J - e s n b (u )} 

2 l-L(« )(__ 1 ) CT (uo) x ( Uo )_ (21) 

This is a 'connected' contribution of order N°. 
This leaves us with the bridging bonds (see Figure |^), 
which may give a nontrivial contribution. However, since 
the bridging bonds occur only for L(uq) > 3 we neglect 
their contribution altogether. 

The subtracted contribution, — i?o(V'a|V ; &)j mav be 
rewritten as 

2 _(jV_i( Uo )) 

-e y^2 1 ~ i(tio) (-i) CT(tio) = 

Uq 

£ 2 p (q+6) -i Aq+Mo; (22) 

a . b 
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where we have used the expanded form ( |I3| ) of x( u o)- 
This is again an energy of order N, which cancels the 
bulk term ( |20| ) up to terms of order N°. After some 
reshuffling (using L(u)=N) we can rewrite the contri- 
bution to the matrix clement as 



(23) 



with 



ej(u) = J (1 — S L(uo)A ) 2 1 - L ^(-iy^ X (u) 



eo(«) 
E{u) 



E{u) - L{u) 



a, 6 uGa+l) 



TV 



(24) 



The first term on the rhs of ( p3|) is the 'bare' on-site en- 
ergy of the triplet. It is always proportional to the over- 
lap integral, so upon switching to the effective Bosons 
this term becomes a fc-independent constant shift. The 
quantity e n (u) may be thought of as describing a 'loss of 
resonance energy'. It is the difference of the two contri- 
butions (20) and ( |2^ ) and occurs because the loop uq is 
fixed, whence the area covered by this loop is not avail- 
able for resonating between different singlet coverings. 
For example we have 



J + eo(l) 



1pa\tpa) 



(25) 



i.e. this quantity is an additive renormalization of the 
energy of formation for a single triplet due to its being 
'embedded into the singlet soup'. Numerical evaluation 
in a cluster shows that this additive correction is quite 
large - for one triplet in a pure singlet background we 
find eo(l) ~ 1.2 J. While this is surprising at first sight 
it should be noted that a similar large value (« 0.8J) 
was previously found in spin ladders [Bj. A fixed triplet 
obviously leads to a quite substantial loss of resonance 
energy. The numerical values of eo{u) and Eq were again 
obtained by cluster calculations (Appendix). Introduc- 
ing eo(it)=eo(w) + Jx( u ) we can write down the third part 
of the Hamilton matrix: 



e< 3 >(fc)=e (l) + £ 1^74*0 

L>2 



(26) 



We can now add up the three contributions, (|16|), (18), 
and ( |26] ) to obtain the total 'connected' Hamilton ma- 
trix e(fc). This is still expressed in terms of the non- 
orthogonal orbitals ft . What remains to be done there- 
fore is to transform the Hamilton matrix to the orthog- 
onal orbitals (|l2). To that end we take matrix elements 
of the form (y, k e(fe)|z/, fe). Introducing the 2x2 trans- 
formation matrix 



(27) 



the transformed Hamiltonian then can be expressed as 
e{k) = T T e{k)T. 

IV. PAIR CREATION AMPLITUDE 

In this section we proceed to a discussion of the pair 
creation amplitude in equation (JoJ) . As discussed in 
section II, by starting from an arbitrary singlet covering 
of the plane and acting with the Hamiltonian along a 
bond connecting two different singlets, we only create a 
state where both singlets are replaced by triplets, see the 
first of Eqs. (|). 

The situation where the two singlets 
in question are 'parallel' to each other, i.e. that the 4 
sites belonging to the 2 singlets form a square with edge 
1, requires special attention. As discussed above we have 
the identity 



S i,i' S j,f\ 



= t 



(28) 



i.e. the state with two parallel triplets can be expressed 
as a linear combination of the two perpendicular com- 
binations of parallel singlet states. In other words: this 
state is already exhausted by the singlet background, and 
consequently must be omitted from our reduced Hilbert 
space. The triplets thus have to obey the additional con- 
straint of never being parallel to each other - one impli- 
cation is that we must set the respective pair creation 
matrix element to zero. 

To compute the numerical value of the matrix element, 
let us consider the action of H on the RVB state (||). We 
have 



H\RVB) = Ne s \RVB) 



j(_iy{m,n) ni 



ma.na/ • 



(29) 



m.n a 



The first term on the rhs originates from processes where 
the Hamiltonian has 'hit' a bond covered by a singlet, 
the second one originates from processes where the ex- 
change has acted along bonds connecting two singlets on 
bonds to and n (the prime on the sum indicates that only 
pairs of bonds connected by a bond are summed over). 
The modulus of the respective matrix element is J/4 (see 
equation (]?])) and there is an extra sign which depends 
on the relative orientation of the bonds to and n. The 
dependence of this sign on the orientation is shown in 
Figure |^. Also, we have approximated the normalization 
factor which is included in the definition of \$> m a,na) by 



1 



1 



n(ma,na) ^fh~{ 



(30) 
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The factor of 1/y/n remains from the definition of (||). 
To obtain the pair creation matrix element we now form 
the overlap between (^) and the state „ri fe JO). 



+ 



o 



+ 



o 6 
o o 
o o 



(a) 



O Cy-vO Ot-tO O 
Oj-— jO O O °J^J° 

O O^-yO ° 

o o o o o o 
(b) 



However, this overlap is rather small: in the case shown 
in Figure (0) it is for example — \/3x(3)/8, and obvi- 
ously this is the most unfavorable case. For other relative 
orientations of the two triplets the overlap can be only 
achieved by a loop of length 8, whence these overlaps 
are oc x(4) <C 1. The non-orthogonalities thus are much 
more benign that those for the spinon-states, whence ne- 
glecting them altogether (as we will do henceforth) is 
probably quite justified. 



EXTRAPOLATION TO FINITE TRIPLET 
DENSITY 



FIG. 8. The sign for all different pairs of bonds 

on which pair creation is possible. The bond m (dashed) is 
kept fixed, the bond n is labeled by (-l) CT ( m >"). 

Defining the 2x2 matrix 



C 



*lxx(&) = — 2cos(2fc x ) — 4cos(fc K ) cos(fc y ), 
t' l yy (k) = — 2cos(2fc y ) — icos(k x ) cos(ky), 

t[ !Xy (k) =4sin(^) S in(|) + 4 S m(^) S in(^) ! (31) 

and bearing in mind that (f fc \ft a ,) — N a ^ a i(k), we find 
for the pair creation matrix: 

A{k) = T T N(k)A(k)N T (k)T. 

This completes the derivation of the single-particle terms 
of the triplet Hamiltonian. Before we proceed, let 
us briefly return to the problem of non-orthogonality. 
Strictly speaking, the state with two triplets is not or- 
thogonal to the singlet background either. The reason is 
that if one draws a loop passing through both triplets, 
the time reversal parities of the two triplets cancel, and 
the state has a nonvanishing overlap with a state where 
the loop is covered only by singlets (see Figure ||) . 



o o 



6 b 



(b) 



FIG. 9. A state with two triplets (a) can have a nonvanish- 
ing overlap with a state consisting exclusively of two singlets 
(b). 



In the preceding sections we have computed the various 
overlap integrals and matrix elements for triplet prop- 
agation, pair creation and interaction. In all of these 
cases the matrix element could be factorized into a con- 
tribution from a 'local' transition between different sin- 
glet/triplet coverings along a single or two neighboring 
loops, and a factor which describes the embedding of 
these active loops into the singlet background. Thereby 
we have always given expressions for these 'embedding 
factors' which are valid in the limit of vanishing triplet 
concentration, i.e. we have computed them as they would 
be for a pure singlet covering of the system. Clearly, this 
is inappropriate for the real system, where quantum fluc- 
tuations have admixed a finite density of triplet Bosons. 
In the following we want to discuss how we have to mod- 
ify our theory to take the effect of a finite triplet con- 
centration into account. It should be stressed from the 
very beginning that this is quite obviously a very complex 
problem and we will be forced to apply some relatively 
crude approximations. 

If we want to derive single-particle like matrix elements 
for finite triplet concentration we should consider overlap 
integrals or matrix elements of H of the type 



JlCt! ,l2Ct 2 , 



I* 



(32) 



i.e. m — 1 triplets stay unchanged, and only a single one 
(which without loss of generality can be taken to be the 
first one) changes its position (but maintains its spin). 
In other words, we should calculate the embedding fac- 
tors for singlet coverings containing a certain number of 
'inert' triplets. Thereby we actually have to make the 
major assumption that the matrix element does not de- 
pend significantly on the positions of these inert triplets 
- otherwise, the very idea of a single particle-like propa- 
gation of the triplets would be invalid. One might then 
expect that approximate values can be obtained by dis- 
tributing the m inert triplets in all possible ways (we 
call the number of possible distributions D(m)) over the 
system, and taking the average of the respective matrix 
elements computed for all D(m) possible distributions. 
In this way, the embedding factors acquire a dependence 
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on the density of triplets. 

The numerical calculation in a finite cluster then pro- 



ceeds as follows: we choose to bonds ix,t2- 



which 



obey the various constraints on the relative positions of 
triplets, and evaluate the ground state norm and energy 
according to 



^ p{a+b) - N U A - 



ct-\-b,iv 



(33) 



The calculation of the various embedding factors then 
proceeds in an entirely analogous fashion i.e. in (13), 
(|H)|), and (M) we replace 



a~\-b,i u 



(34) 



a, b 



a , b 



In this way we obtain all embedding factors for fixed dis- 
tribution of inert triplets, and the value for triplet con- 
centration m/N is obtained by averaging over all allowed 
distributions of the bonds ii,«2 •• -im- ln practice this 
calculation requires quite a substantial numerical effort 
so this was performed only for the 4x4 cluster. 



smaller and smaller. Since longer loops may cause prob- 
lems in the small clusters, a suppression of these may 
therefore even reduce finite-size effects. Moreover, the 
suppression of longer loops is actually beneficial for our 
entire approach: it tends to eliminate the problems with 
the nonorthogonalities in the singlet soup and makes the 
truncation of the length L(uq) a better approximation. 
As an illustration Figure [l0| shows the dependence of 
eo(l) 7 i-e. the additive renormalization of the energy of 
formation of a triplet (see (p5f), evaluated in the 4x4 clus- 
ter. The concentration dependence is as expected: for 
low triplet concentration e (l) is large because an added 
triplets blocks many long loops along which resonance 
could have occurred. As the triplet concentration gets 
higher, these longer loops a blocked anyway, so adding a 
further triplet does not increase the energy too much any 
more. In the high-density limit the additive renormaliza- 
tion approaches zero, as expected. The Figure shows, 
however, that the concentration dependence is quite sig- 
nificant, i.e. this effect should not be neglected. 
In addition to reducing the importance of longer loops, 
for a finite density of triplets we have to take care of the 
excluded volume constraint which the bond Bosons have 
to obey. Placing a triplet on one given bond m blocks a 
total of 9 other bonds, on which no more triplets can be 



placed (see Figure 11) 




FIG. 10. Dependence of the various parameters on the 
triplet density 8 as obtained from numerical calculations on 
the 4x4 cluster. 

One might expect, however, that the finite size effects 
are in fact smaller for the more relevant higher triplet 
concentrations: the main effect of the fixed triplets ob- 
viously is to reduce the importance of long loops. Since 
the bonds occupied by the static triplets must be iden- 
tical in the bra and ket state, loops of length L(u) > 2 
which pass through these bonds are impossible. When 
the density of triplets gets appreciable, the probability 
to find enough space for forming longer loops becomes 




FIG. 11. Placing a triplet on the central (dashed) bond 
blocks the 9 indicated bonds, on which no more triplet can 
be placed. 

Of these, 7 bonds are blocked because they share a site 
with m, the remaining 2 are blocked because they are 
'parallel' to m and the state with two parallel triplets 
is actually a linear combination of singlets (see (§8h). 
In order to take care of this blocking effect, we resort 
to a Gutzwiller-type approximation. It has been shown 
by Ogawa et al. [|TJ that the essence of the Gutzwiller 
approximation is the neglect of the difference in phase 
between states where the particles in question are dis- 
tributed in different ways over the lattice. With this 
approximation, any real space distribution of m Bosons 
contributes the same number (which without loss of gen- 
erality can be taken to be 1) to the norm of any state 
with m Bosons. The total norm of any such state then 
becomes simply the number of possible ways to distribute 
to Bosons over the plane, i.e. D(m). In the present case, 
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the situation is somewhat more complicated, because we 
still have to take into account the fact that due to the 
loop problem the norm of a state with to fixed triplets 
is not 1 (unlike the case if the triplets were just ordinary 
Bosons). We thus approximate the normalization factor 
of any state with to triplets by 

n ( m ) = /r } ; _ > ( 35 ) 
y/D[m)n m 

where n m is the value obtained by averaging (|33| ) over all 
allowed triplet configurations ii,%2 ■ ■ - im- Clearly, this 
replacement will lead to a renormalization of the various 
embedding factors. To evaluate these, we usually have to 
consider a certain area A, covered by one or two loops, in 
which the actual overlap, hopping or pair creation pro- 
cess takes place. We then determine (again by simulation 
on the 4x4 cluster) the number of ways to distribute m 
triplets over the exterior of the area A. Thereby we re- 
quest that putting a triplet anywhere within A always 
gives an allowed configuration of to + 1 triplets. We call 
the number of allowed triplet configurations D(m,A). 
Defining the various areas Ai as 



To summarize this section, Figure |l0| shows the values of 
the different parameters as estimated by the procedure 
outlined above, for all possible triplet concentrations in 
the 4x4 cluster. Let us stress again that the approxima- 
tions leading to the parameters in Figure [l^ are proba- 
bly rather crude - one may expect, however that we get 
roughly correct orders of magnitude and that the ratios 
of the different parameters come out approximately cor- 
rect. 



VI. SPIN DYNAMICS 

Let us first briefly summarize the discussion so far: 
We have shown that the elementary excitations of the 
singlet soup correspond to excited dimers, which must 
be modeled by bond-Bosons. Combining the matrix ele- 
ments computed in the two preceding sections we obtain 
a Hamiltonian of the form 
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FIG. 12. Areas on which no triplets are allowed in some 
processes. 

in figure [l^, we then have to renormalize the various pa- 
rameters as follows: 

ni 



X(2) - X(2) 
co(2) -> e (2) 



y/D(m,Ai)' 
D(m,A 2 ) 



D(m,Ai) ' 
D(m,A 2 ) 



T) — > T) 



D(m,Ai)' 
D{m,Ai) + D(m,A 5 ) 
2D(m,Ai) ' 
D(m,A & ) +D(m,A 7 ) 



2D(m,Ai) 

For the pair creation amplitude we replace 



C D{m,Ai) 



D(m,A 5 ) 



(36) 



(37) 



Thereby the matrix elements A(fc) and e(fe) are func- 
tions of the triplet density. For given triplet density Si 
these parameters can be computed and the Hamiltonian, 
which after the Gutzwiller-type renormalization of the 
matrix elements we take to be a free-Boson Hamiltonian, 
is solved by Bogoliubov transformation. Combining the 
two different r-operators into a 2-vector T the ansatz 
reads 

T k =u k T k +v k Tl k , (38) 
where the real 2x2 matrices u and v have to fulfill 



u.uT — VuvT = 1, 



M k u k v k u k 

u k v k ~ v k u k = 0- 



The inverse transformation of (B8[) is therefore 



Tk — u k^k v h^-k' 



(39) 



(40) 



The Hamiltonian can be transformed to free particle 
form: 



+ 3^2 Tr[v h e(k)vl - v k A(fc)uJ - u k A(k)vl] , (41) 



provided the transformation matrices obey 

e(k)ul-2A(k)vl = uluj(k), 
2A(k)ul - e{k)vl = vluj(k). 



(42) 
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Here we have introduced the 2x2 matrix uj{k) 
diag(uJi(k),LU2(k). The triplet density is 

3 



Tr[u k u k f{uj) 



/H)], (43) 



where / is the Bose function. The requirement Si = Sf 
then provides a self-consistency condition for the density. 
Numerical evaluation shows that there is a minimum 
triplet concentration <5o- For 6 < Sq, the parameters 
are such that in a certain area around k — (jr, tt) the 



eigenvalue problem (42) does not have real eigenvalues. 
For 8q the minimum of the triplet dispersion, which oc- 
curs at (tt, 7r) is precisely zero. With the coarse mesh 
of concentrations possible in the 16-site cluster, we find 
0.25 < Sq < 0.375. We therefore linearize all matrix ele- 
ments x, i.e. x(S) = a x + b x {5 — 0.25), using the values 
at 5 = 0.25 and S — 0.375 to determine a x and b x . This 
gives S = 0.335. 

With increasing temperature, the self-consistent value of 
S increases slowly, approaching So for T — ► 0. Figure 
|l3| shows the dispersion relation for the respective self- 
consistently determined triplet density at different tem- 
peratures. Since we have two degrees of freedom/site 
(the triplet in x and y direction) we obtain two bands. 
While one of these bands is practically dispersionless, the 
dispersive one resembles results obtained for other spin 











f?J=2 












JSJ=20 \ 







(0,0) 



(IF,*) (ff.0) 



(0,0) 



FIG. 13. Dispersion of the two eigenvalues u)i, for two dif- 
ferent temperatures. 

liquids, such as the two- leg ladder |8||| or the bi-layer 
Heisenberg antiferromagnet H] : the dispersion starts at 
relatively high energy at k = (0, 0), takes a shallow max- 
imum near the antiferromagnetic zone boundary and, as 
stated above, takes a more or less pronounced minimum 
at Q = (w,ir). For low temperatures the bandwidth is 
w 2 J, which is quite close to the value for antiferromag- 



netic spin waves. The case of zero temperature requires 
special attention: it is obvious from Figure [l^ that the 
gap at Q approaches zero for T — > 0. Assuming that the 
gap vanishes at zero temperature, we may then assume 
that this momentum becomes macroscopically occupied 
by a triplet density Sq. To determine this condensed frac- 
tion we note that the limiting density Sq is defined such 
that the gap at Q is exactly zero with the parameters 
calculated for this density. Then, for Si = S we evaluate 
the density Sf of uncondensed triplets (i.e. we exclude 
k = Q from the sum in (E3)) and finally determine the 
condensate density from 



So 



(44) 



This is entirely analogous to the treatment of Hirsch and 
Tang |2^] of the condensation of Schwinger-Bosons in the 
mean-field theory of Arovas and Auerbach |2j|. As will 
be seen in a moment, just as for the Schwinger Bosons 
the condensation of the triplets corresponds to antiferro- 
magnetic ordering. The ground state energy then is 



E tot = E (S) 



3j]Tr[»nt)fw(t)] 

k 



(45) 



where Eq(S) is the energy of the singlet background for 
triplet density S. This is calculated from (^). We ob- 
tain a value of — 0.3426J/bond, which is quite close to 
the true ground state energy of the 2D Heisenberg anti- 
ferromagnet |2^ 1 . An interesting figure is the lowering of 
the energy as compared to the original RVB- vacuum (^) . 
In the 4x4 cluster (which has been used to compute 
all matrix elements) the expectation value of the pure 
nearest-neighbor RVB state is — 0.334318J/bond [||, so 
that the admixture of the triplets lowers the energy only 
by a tiny 0.008 J/bond. In view of the strong approxima- 
tions we were forced to make this result probably has lit- 
tle quantitative significance - it shows quite clearly, how- 
ever, that the energy of the RVB-vacuum is lowered only 
by a very small energy by the triplet fluctuations. This is 
what must come out because in the thermodynamic limit 
the energy of the RVB state is —0.302 J/bond ]lt^| , com- 
pared to the exact ground state energy of approximately 
0.334 J/bond @. 

We proceed to a discussion of the spin correlation func- 
tion. Since the bond-Bosons are actually objects which 
extend over more than one unit-cell, they have something 
like a structure factor which depends on the type of op- 
erator by which they are probed. Let us consider a single 
bond, and introduce even and odd combinations of 

spin operators on this bond: 



Si 



Sf ± tf? 



Due to their opposite parity under exchange of i and j 
these two operators have an entirely different effect on 
the 4 possible spin states of the bond: S+ annihilates 
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-it* 



the singlet and t\ but converts tj, — > zi y and ft 
whence ST^t* x t. S z _ on the other hand converts the 
singlet into t\ and vice versa, but annihilates both t\ and 
Si, whence S f ^=r^ + r. If we still restrict ourselves to 
only a single bond n we have to 'translate' 



( — (n + t„)], (46) 

where <7„ is the unit vector along the bond n and R n 
denotes its center of gravity. Let us next discuss how 
these matrix elements are modified due to the embedding 
into the singlet background. We act with the operator 
S(q) = (1/V2N) Ej e iq - R JSj on some given state. The 
first term simply changes the z-spin and momentum of a 
triplet. It does not affect the number of triplets, so we 
need no further renormalizations. For the second term, 
let us start out from a state with one triplet, \j, a), and 
assume that the triplet is annihilated, i.e. converted into 
a singlet. This leaves us with the state 



1 S 



N-l 



The matrix element with the pure RVB state then is 



(RVB\^) = 



1 



1 



because it is easy to see that a state with one fixed singlet 
has an overlap of n/A with the RVB state. If we want 
to extrapolate this to finite triplet density we again have 
to use the values of n and n\ computed with m fixed 
triplets, and renormalize the matrix element by 



! D(m,Ai 
D{m) 



We thus find 



-j^Bm(^)(Tt ia + r, ia )]. (47) 

By using the inverse transformation (^) this can now 
be converted to the eigenvectors T. The dynamical spin 
correlation function of the spin liquid thus consists of 
two components: First, there is a two-particle contin- 
uum, which dominates for momenta around (0,0). Sec- 
ond, there is a single-particle like contribution, which 
dominates the cross section for momentum transfers near 
(n, 7r), and hence should be identified with the excitations 
seen in neutron scattering around this momentum. The 
situation is the same for ladders, where the single-particle 



spectrum is observable in the channel with momentum 
transfer perpendicular to the ladder, k± = 7T, and the 
two-particle continuum for k±_ = M. 
Next, we discuss the relationship between condensation 
of triplets and antiferromagnetic ordering. The opera- 
tor of staggered magnetization (which is a vector) can be 
written as 




T 0,«) 



(48) 



Here we have used the fact that for the high-symmetry 
momentum Q the overlap and Hamilton matrix are triv- 
ially diagonalized by the symmetric and antisymmet- 
ric combinations of the bonds in x and ^/-direction , 

X Q,± = TS^Q,* ± f Q,y)- Then ' introducing the 3D 
unit vector fi, we can construct the coherent state 



XVNCl- 



'Q>+\RVB). 



This state corresponds to a condensate of the Tq + 
Bosons, whose density is given by Sq = A 2 . The stag- 
gered magnetization per site is 



K 5 Q n(l + xW) 
2 V ni 



n. 



Inserting the value for the condensate density 5q ob- 
tained from (Q) we obtain the value m s = 0.25. Most 
current estimates for the 2D Heisenberg antiferromagnet 
are around m s = 0.35 |24[ . 

As mentioned above, the identification of antiferromag- 
netic ordering and condensation of some kind of effec- 
tive Bosons is quite reminiscent of the treatment of 
Hirsch and Tang [p2| in the framework of Schwinger- 
Boson mean-field theory In the present theory the 
direction of the staggered magnetization is given by the 
unit vector ft, which can be chosen arbitrarily. Conden- 
sation of the triplet Bosons determines the total density 
of triplets, but not their distribution over the three spin 
species. The 'ground state' thus is actually an entire 
manifold of states which can be transformed into one an- 
other via SO(3) rotations of the vector ft. One might 
thus conjecture the existence of low-energy states where 
the direction of ft changes slowly in real space. These 
states, which corresponds to a slow fluctuation of the an- 
tiferromagnetic order parameter may be describable in 
terms of a nonlinear a-model - wc defer a detailed dis- 
cussion of this issue to a separate paper. 
To conclude this section, we discuss the relationship with 
Zhang's SO (5) symmetric theory of cuprate supercon- 
ductivity [lq ]. The preceding discussion has shown that 
an antiferromagnetic state can be viewed as a spin liq- 
uid where the triplet-like bond Bosons have condensed 
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into the state which has momentum (ir, tt) and s-like 
symmetry under point group operations. The antiferro- 
magnetic order parameter, a real 3-vector, then is the 
vector of condensation amplitudes of the three bond- 
Boson species. This interpretation of the antiferromag- 
netic state in fact is a key ingredient for a microscopic 
interpretation of Zhang's SO(5) theory of cuprate super- 
conductor. Namely the 7r-operator, which acts as the 
'ladder operator in charge direction' in Zhang's repre- 
sentation of the SO(5) angular momentum algebra Jl5| , 
precisely converts an s-like combination of bond triplets 
with momentum (tt, tt) into a nearest neighbor e?-wave 
hole pair. This can be seen by writing the 7r-operator in 
real space: 

TT* = ^V Q ' fll [ ( Cijd+i^ + c id c i+x ^ ) 
i 

-( Ct,lCi+y,l + Cid c i+y,1 ) ]■ ( 49 ) 

Introducing a further bond Boson h\ which stands for 
the hole pair, this could be written as 

where the momentum K = (0, 0). Acting with (f2 • ir) pN 
onto an undoped state with the triplet condensate den- 
sities pfi will therefore convert this state into a conden- 
sate of d-like nearest neighbor hole-pairs with momentum 
(0,0). Viewed in this way, the connection between the 
antiferromagnetic and the superconducting state appears 
quite natural. The only drawback would be that the hole 
pairs are strictly nearest-neighbor pairs, i.e. the present 
version of the 7r-operator neglects charge fluctuations. 



VII. DISCUSSION 

In summary, we have discussed the excitation spec- 
trum of a completely disordered and homogeneous spin 
system. Thereby we used an approach which might be 
viewed as a generalization of spin wave theory: whereas 
spin wave theory describes fluctuations around a Neel or- 
dered 'vacuum', we have instead used the nearest neigh- 
bor RVB state as basis for self-consistently constructing 
the excitation spectrum. This general idea of treating 
a completely disordered state is probably more widely 
applicable, for example to treat the 'orbital liquids' pro- 
posed recently Eq] for manganates. 
As a first key result, we then found that the elemen- 
tary excitations of the singlet soup-like vacuum are not 
Fermionic 'spinons', but rather Bosonic excitations Tj, 
which can be viewed as excited dimers propagating 
through the system while constantly resonating between 
x and y direction of the dimer. Following a similar pro- 
cedure as in the Zhang-Rice derivation of the t-J model, 



we could write down a Hamiltonian for these Bosonic ex- 
citations. The Hamiltonian contains an energy of forma- 
tion for the triplets, a term describing their propagation 
and a term describing pair creation and annihilation. All 
quantities in the derivation are quite well-defined, and 
can in principle be obtained by numerical techniques. In 
the present work we have actually attempted an 'ab ini- 
tio calculation' of the various parameters in the effective 
Hamiltonian - a promising alternative to this somewhat 
clumsy approach might be a semi-empirical approach, 
where the parameters in the Hamiltonian are adjusted to 
match e.g. experimental data. 

The obtained ground state is an exact spin singlet, trans- 
lationally invariant and isotropic - in short, it has pre- 
cisely the symmetry properties expected for a spin liq- 
uid. The elementary excitations is a triplet mode, which 
reaches its lowest energy at Q=(tt,tt). Such a triplet 
mode is the most natural generalization of an antiferro- 
magnetic spin-wave to the spin-liquid state. Condensa- 
tion of these bond Bosons into the state with momentum 
(tt, tt) then corresponds to antiferromagnetic ordering of 
the system. The breaking of rotational symmetry in spin 
space is due to fixing (different) condensation amplitudes 
for the three components of the triplet mode. The latter 
correspond to the three possible components of the anti- 
ferromagnetic order parameter. 

In the present work we have chosen the completely 
isotropic nearest neighbor RVB-state as the starting 
point for constructing the triplet-Hamiltonian. This is 
appropriate for an truly isotropic system, which probably 
is realized at finite doping and/or sufficiently high tem- 
perature. Using a translationally invariant state, how- 
ever, is not mandatory. For example by redefining the 
operator S as 

s = £((i + e)4^* + (i-e)4, i+fi )> ( 5 °) 

i 

with e > we can obviously generate a 'singlet soup' 
with a preference for singlet orientation in x-direction. 
It is tempting to speculate that this may be appropri- 
ate for orthorhombic Lai_ a ;Sr x Cu04, where at x = 0.125 
columnar order is known to exist. Analogous calculation 
as above, but with a finite e may thus be appropriate to 
discuss the excitation spectrum of this material. 
Finally we note that the present scenario for the excita- 
tion spectrum of the spin-liquid provides a corner stone 
for a microscopic interpretation of Zhang's SO (5) sym- 
metric theory of cuprate superconductors pa |. As dis- 
cussed above, an antiferromagnetic state may be thought 
of as being generated by condensing triplet-like spin ex- 
citations of the RVB spin-liquid into the state with mo- 
mentum (jr, tt). The antiferromagnetic order parameter 
corresponds to the vector of condensation amplitudes of 
the three triplet species. Then, SO (5) symmetry states 
that such a triplet-excitation with momentum (tt, tt) is 
'dynamically equivalent' to a d x i_ y i hole pair p5|. In 
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fact the 7r-operator |15j, which plays the role of a ladder 
operator for 'rotations in charge-direction' in the SO (5) 
theory |fL5), precisely replaces a bond-triplet with mo- 
mentum (7r, 7r) by a hole pair with momentum (0, 0). The 
fact that the 7r-operator is an approximate eigenoperator 
of the Hamiltonian, [H, n] = lo^t:, then implies that the 
triplet and the hole pair are dynamically indistinguish- 
able and differ only by their energy of formation. The 
7r-operator thus would convert a condensate of triplets 
into a condensate of hole pairs, and, by the dynamical 
equivalence of these two objects, thus converts the anti- 
ferromagnetic ground state at half-filling into a supercon- 
ducting state at finite-doping one. The somewhat 'toy- 
model-type' discussion of Ref. |^5| thus may very well be 
transferable almost literally to the fully planar t-J model. 
In the present work we have restricted ourselves to a pure 
spin system. One my expect, however, that all considera- 
tions go through also for the doped system, with the sole 
difference that we get additional renormalizations of the 
matrix elements in the effective Hamiltonian due to the 
fact that the doped holes reduce the volume available for 
pair generation and propagation of triplets. Moreover, 
we will need terms which describe the coupling of the 
triplet branch to the mobile holes. In the simpler case 
of hole motion in a spin ladder this program has in fact 
already been carried out ||^7| , leading to quite satisfac- 
tory agreement with numerical results. For the planar 
case we defer this to future work. 

I would like to thank W. Hanke, O. P. Sushkov and Shou- 
Cheng Zhang for instructive discussions and helpful com- 
ments. 



Moreover, the data show that the x(L) an( i decrease 
quite rapidly with L. truncation of the series after L = 2 
thus seems to be quite a reasonable approximation as 
well. 





M=2 


M=4 


L 


m 


ML) 


X(L) 


ML) 


1 


0.10953 


0.13844 


0.121483 


0.15142 


2 


0.03602 


0.01675 


0.043671 


0.01998 


3 


0.00780 


-0.00197 


0.012879 


-0.00001 



TABLE I. The values of x(L) and e (L) for the two differ- 
ent L. 



VIII. APPENDIX 



In this appendix we discuss the numerical work needed 
to obtain the values of x( u )-> e o( u ) an d t ne rj. To that 
end all possible dimer coverings of an M x M cluster with 
periodic boundary conditions are generated on the com- 
puter the sums in (13) are evaluated numerically. Re- 
strictions on computer memory and CPU time do not 
allow to use M > 6, so that in practice only M = 4, 6 
are possible. Tu study the size dependence, we consider 
the quantities = — x(£); these give the ratio of the 
norm of an RVB state covering the exterior of a loop of 
length L to the norm of an RVB state covering the en- 
tire cluster. Similarly, we define eo(L) = ^eo(L), which 
gives the gain or loss in energy due to a fixed loop of 
length L. The values for the different x an d £o are given 
in Table [j] and actually show already a quite remarkable 
independence of systems size. Moreover, the expecta- 
tion value of the energy/bond, Eq/N, is —0.334318 for 
M = 4, —0.313763 for M — 6; the estimate for L = oo is 
—0.302 |l6). All in all the data suggest that already the 
4x4 cluster gives reasonably accurate estimates for the 
different parameters. 
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one could also construct a basis of 'spinon' states and 
re- interpreting 
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